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1.  INTRODUCTION 

It  is  well-known  that,  classical  elasticity  solution  of  the  crack  and  dislocation 
interaction  possess  singularities  both  at  crack  tips  and  at  the  point  of  application  of  the 
dislocation  (cf.  [1]).  The  non-physical  nature  of  this  solution  has  been  debated  exten¬ 
sively  in  the  literature  until  few  years  ago,  when  the  solution  of  the  same  problem,  based 
on  nonlocal  elasticity,  appeared  (see  [2]).  According  to  the  nonlocal  elasticity  the  stress 
field  is  finite  at  crack  tips  [3,4]  and  at  the  dislocation  core  [5,6]  but  possess  maxima  near 
these  points.  Thus,  a  plastic  region  constructed  by  means  of  an  appropriate  distributions 
of  dislocations,  will  not  contain  any  stress  singularity  [7]. 

According  to  physics  of  solids,  the  initiation  of  fracture  is  attributed  to  the  max¬ 
imum  stress  field  exceeding  the  cohesive  stress  that  holds  atomic  bonds  together.  Conse¬ 
quently,  nonlocal  results  permit  the  introduction  of  the  maximum  stress  hypothesis  as  a 
fracture  criterion  [4,8],  and  this  is  shown  to  give  excellent  agreements  with  atomic 
theories  and  experiments  [9,10]. 

Motivated  with  these  results,  here  we  pursue  to  discuss,  a  primitive  problem  by 
means  of  which  one  can  construct  the  solution  of  the  small  scale  yielding  in  the  forms  of 
two  inclined  straight  lines  emanating  from  each  crack  tip.  This  model  is  an  accepted  one, 
both  on  the  basis  of  experiment  and  classical  treatments  [11,12].  This  model  can  be  used 
to  construct  solutions  of  plastic  yielding  based  on  distributions  of  dislocations  on  areas 
near  the  crack  tips.  Such  a  model  is  much  more  realistic,  and  it  is  in  accordance  with 
experimental  observations. 

The  solution  of  Mode  ID  crack,  with  a  line  distribution  of  screw  dislocations  along 
the  crack  line,  was  presented  in  a  previous  work  [13].  The  present  paper  represents  an 
extension  of  the  same  problem  to  the  Mode  I  crack,  with  inclined  edge  dislocations. 
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2.  BASIC  EQUATIONS 

Linear  theory  of  nonlocal  elasticity,  for  homogeneous  and  isotropic  solids,  (vanish¬ 
ing  body  force  and  inertia)  is  governed  by  Cauchy’s  equations  of  motion 

tu.k  =  0  (1) 

and  the  stress  constitutive  equations,  (cf.  [2,9]), 

tu  =  Jctflx'-x \)au(x')dv(x')  (2) 

v 

where  is  given  by 

o«  =  Xurtr^ki  +  +«/,*)  (3) 

Here  tu(x)  is  the  real  stress  tensor  at  the  reference  point  x.  It  is  influenced  by  the  strains 
at  all  points,  x',  of  the  body,  with  volume  V  enclosed  within  surface  BV.  X.  and  p  are  the 
usual  Lame'  constants  and  uk(x')  is  the  displacement  vector  at  x'eV.  The  influence  func¬ 
tion  a<ljc'-x  I)  depends  on  the  distance  lx' -x  I  of  x'  from  x.  Various  forms  of  this 
function  was  given  by  Eringen  [14].  Here  we  use 

a(lx  1)  =  (2ite2)' 1  K ^  /  e)  (4) 

which  is  appropriate  to  the  plane  strain  problem.  Here  Kq  is  the  modified  Bessel’s  func¬ 
tion  of  the  first  kind,  and  e  is  an  internal  characteristic  length  appropriate  to  material; 
e.g.,  for  perfect  crystals  it  can  be  taken  as  the  lattice  parameter.  For  other  kernels  in  1, 2 
and  3  dimensions  see  [14], 
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It  is  of  interest  to  note  that  (4)  is  the  Green’s  function  of  the  linear  operator 
L  = 1  -e2V2,  i.e., 


(1  -e2V2) a(lx'-x  I)  =  5(lx'-xl) 

where  V2  is  the  two-dimensional  Laplacian  operator  and  Ux)  is  the  Dirac  delta  measure. 
This  feature  of  the  function  K  o  allows  us  to  invert  the  constitutive  equations  (2)  of  nonlo¬ 
cal  elasticity: 


(1  -  e2V2)  tu  =  ou  (5) 

With  this  apparatus  at  hand,  Eringen  [15]  gave  the  solution  of  the  problem  of  the  stress 
field  due  to  edge  dislocations  which  is  used  in  the  following  analysis. 
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3.  FORMULATION  OF  THE  PROBLEM 

The  main  purpose  of  this  study  is  to  describe  the  stress  field  in  an  infinite  elastic 
plane,  weakened  by  a  line  crack  and  has  four  single  dislocations,  which  are  located  sym¬ 
metrically  with  respect  to  the  crack,  Fig.  1.  The  crack  is  subject  to  a  constant  tensile 
stress  t 0  and  the  Burger’s  vectors  of  the  four  edge  dislocations  are: 


-  b  (cos  a,  sin  a,  0) ,  -b(-  cos  a,  sin  ot,  0)  (6) 

-  b  (cos  a,  -  sin  a,  0) ,  -  b  (-  cos  a,  -  sin  a,  0) 

The  stress  field  is  determined  by  the  superposition  of  three  primitive  stress  fields  in 
the  medium: 

(i)  Stress  field  due  to  the  constant  traction  on  the  crack  surface  (r  ®) 

(ii)  Stress  field  due  to  the  dislocations  (r£) 

(iii)  Stress  field  due  to  the  interaction  of  crack  with  dislocation  (tru) 

The  solution  of  problem  (ii)  was  given  by  Eringen  [15].  The  stress  field  due  to  an 
edge  dislocation  with  Burger’s  vector  b(  1,  0,  0),  and  located  at  the  origin  of  the  coordi¬ 
nates,  can  be  expressed  as  follows: 

(x,y)  =  -  k/"3/  i  (p)  sin  0  +  /3  (p)  sin  3  0;  ,  (7) 

$  (x,y)  =  -  K{f  i  (p)  sin  0  -  /3  (p)  sin  3  0;  , 

( x,y )  =  -  K{f  i  (p)  cos  0  +  /3  (p)  cos  0; 


where 
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v  -  Uo  n  -  Z. 

4k(1~  v)e  •  p  e 


! 

I 

i 

(8)  | 


p  oo 

/l  =  /3(p)  =  *3(p)J/3(pW  +  /3(p)  J  ATsCpO  rfp' 

p  o  p 

l 

Here  r  =(x2  +  y2)^  and  0  are  polar  coordinates,  Im  and  Km  stand  for  the  modified 
Bessel’s  functions  of  first  and  second  kinds  of  mih  order. 

The  problems  (i)  and  (iii)  can  be  combined  and  solved  together.  To  this  end  we 
need  to  describe  the  load  on  crack  surface  imposed  by  dislocations.  Considering  the 
geometry  of  the  problem  (Fig.  1)  this  load  can  be  expressed  as  follows: 

o'  =  t„  sin2  a  +  cos2  a  +  2t^  cos  a  •  sin  a  (9) 


where 


=  2[r£(pi,  6,)  +  tg{  p3, 63)] 


(10) 


tyy  —  2[tg(P\,  9l)  +  fyy(p3,  93)] 


=  2u5(pi,ei)  +  ^(P3,e3)] 


Pi  =  [(c  - x)2  +  a2  +  2a (c  -x)cos  a]7 ,  0i  =  w  +  Vi  -  a,  Vi  =  arc  sin 


P3  =  ((c  +x)2  +  a2  +  2a(c  +  x)cosa]7,  03  =  x  +  y3-a,  =  arc  sin 


Pi 


P3 


sin  a 


sin  a 


(ID 


To  solve  a  crack  problem  which  is  loaded  arbitrarily  we  follow  Eringen  [16].  The 
stress  and  displacement  field  for  plane  strain  problems  can  be  expressed  as  follows: 


u(xy) 


v{x,y) 


^2n  -»  4y2 


tyy  =  — L  J  %M^y)e-i^X  dl, 


(12) 


4 - 5-  f  £-N 

^  V 


<2^-4 


jV^M^.y)  +  (i  -v)  -  e2 


v4-V+(2v-l)  <gv+(\-v)ffi 


*-'5*  d% 


1 

■25T 


^]^{<v-2)5^  +  (l-v) 


dV  _  p2 


(y-lV^-  +  (2v-3**  +  (1-v) 


1) 


e"*^  d4 


where 


•p(ty)  =  [^(O  +  y  +C©e*V^>  y£0 


(13) 


The  unknown  function  A  (£),  D(%)  and  C(£)  will  be  determined  by  the  boundary  condi¬ 
tions  for  specific  problem  under  consideration. 

For  a  crack  located  at  y  =  0  ijc  I  <  c  and  loaded  symmetrically  with  normal  trac¬ 
tions,  i.e. 


tyy  =  -OoOO  =  <y(jc)-o°,  (Oo(x)  =  cr0(—  x ))  (14) 

we  have  the  boundary  conditions 

K 


txy  (.X,y)  =  0,  y  -  0 


(15) 
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and 

v(x,y)  =  0 ,  y  =  0  \x  \  >  c  (16) 

The  boundary  condition  (15)  yields 

D£)  =  1^1  A®  +  C<&  (17) 

Moreover,  to  ensure  that  ta  is  finite  at  the  tip  of  the  crack  we  set 

C<£)  =  ^2+e-2)"^  K^  +  e-2)^  -  2^r'A©  (18) 

With  the  help  of  (17)  and  (18),  the  boundary  conditions  (14)  and  (16)  read 

VFi  £B(£)[1  +  AT(eS)]cos(x£)rf£  =  -Oo(x)  lx  I  < c  (19) 

0 

00 

Js(^)cos(x^)d^  =  0  lxi>c  (20) 

0 

Here 

"te2+  l  [1  - v  +  (e£)2]  +  (e4) [- l+v-(e4)2]  „ 

B  <g)  =  ^ - —= - 1 - A  <g)  (21 ) 

V(e4)2+i-2(e4) 

*(«)  =  *i(«)[l  +  K  2(« )]  +  2(«)  (22) 
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Ki(u)  = 


Vtt2+l(VU2+l-2tt) 


,  K*U)  = 


Vu2+  1  (V-M2)  +  lifl-V+li2)  -  2m 
■^U2+  1  (1  -v  +  U2)  +  u(-  Y  +V  -  K2) 


(23) 


After  the  dual  integral  equation,  given  by  (19)  and  (20),  are  solved,  (i.e.,  B(£)  is 
described  so  that  the  equations  (19)  and  (20)  are  satisfied)  the  stress  and  the  displacement 
fields  are  found  from  (12)  taking  into  account  the  relations  (17)  and  (18). 


-10- 


4.  SOLUTION  OF  THE  PROBLEM 

It  is  convenient  to  introduce  non-dimensional  quantities 

*-f.  n  =  3  =  f  (24) 

The  dual  integral  equation  (19)  and  (20)  may  now  be  expressed  as 

oo 

I qS(q)[l+K(3q)]  cos  ,zq)dq  =  -S(z),  lz  I  <  1  .  (25) 

o 

oo 

J  B(q)  cos  (zq)  dq  =  0  lz  I  >  1  .  (26) 

o 

where 

S(z)  =  c2(-J)^ao(cz)  (27) 

The  solution  of  this  dual  integral  equation  is  not  known.  For  an  approximate  solution,  we 
recall  the  expression  [17, 6.693.2] 

oo 

J  T|" 1  /2*  +  i0l)  cos  (zT|)  dq  =  0 ,  (*=  0,1,2,  •••).  lz  I  >  1  (28) 

o 

Accordingly,  if  B  (q)  is  chosen  as 

s 

£*0l)  =  I  bkTr'Ju  +  l  (q) 

*  =o 


(29) 
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then  the  displacement  boundary  condition  (26)  will  be  satisfied.  To  determine  the  unk¬ 
nown  coefficients  bk,  we  introduce  (29)  into  (25) 


N 

SN(z)  =  £  bk  j72k  +  i(Tl)[l+Ar(Pti)]-:s(2ti)dTl  (30) 

k  =0  0 


and  impose  that 


r(z)  =  [S(z)-SN(z)]2  (31) 

be  a  minimum.  To  provide  economy  in  computer  calculations,  we  take  into  account  the 
limit  behavior  of  the  kernel  K  (pp)  in  (25) 

lim  (pn)  ->  (32) 

p  -*o  1-v 

Computer  experiments  indicate  that  for  P<10-2  the  kernel  AT  (Prj)  behaves  like  (32) 
except  near  the  crack  tips.  With  this  value  of  K  (Pp)  a  great  deal  of  simplifications  is 
achieved,  since 


Sk(z)  =  J  /  at  +  i  (H)  cor  zr\)dr\  =  cos  [(2*  +  1)  arc  sin  r]  ,  z<1>  jfe=  0,1,2,  •••  (33) 

0  (l-Z2)7 

see  [17,  6.671.2].  For  integer  values  of  k  these  functions  can  be  expressed  in  the  form  of 
polynomials 


*  +  i 

Ski?)  =  Z  AyX20-,), 


/-I 


(34) 
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where 


i4 1  =  1  ,  Aj  +  i  =  - 


(35) 


Since  r(z)  £  0  for  every  z ,  r(z)  is  minimized  in  the  following  sense 


J  r(z)  dz  -  Az  £  r(zi ),  z,  =  Az  (/  -  1) ,  A z  =  ^  ]  (36) 

o  <  =1 
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5.  RESULTS  AND  DISCUSSION 

By  means  of  the  minimization  process  described  in  Section  4,  the  unknown 
coefficients  B  (rj)  is  determined.  Afterwards,  the  stress  field  along  the  crack  line  is  calcu¬ 
lated  by  use  of  Eq.  (18)  and  (12).  The  crack  is  loaded  by  a  constant  tensile  load  t°  and 
four  dislocations  symmetrically  located  along  two  straight  slip  lines  emanating  from  each 
crack  tip.  Each  slip  line  makes  an  angle  a  with  the  crack  line  as  shown  in  Fig.  1.,  and 
one  edge  dislocation  is  located  along  each  line  at  a  distance  a  from  the  crack  tip.  The 
crack  length  is  2c  and  the  nonlocality  parameter  is  e. 

Calculations  have  been  performed  for  aluminum  with  properties: 

fi  =  2.51  x  10u  cgs,  v  =  0.3 

b  =  4.05  x  10" 8  cm 

In  order  to  see  the  effect  of  the  location  of  dislocation  in  one  set  of  calculations 
(depicted  by  Figures  2-5)  the  distance  a  is  varied  while  a  is  fixed,  and  in  the  other  set 
(Figures  6  and  7)  a  is  fixed  and  a  is  varied;  i.e. 

(i)  a  =  67°;  a  =  10~4 ,  10~5,  5xl0~5,  5xl0-4cm 

(ii)  a  =  60°,  75°;  a  =  10" 4  cm. 

In  calculating  B/v(t|),  given  by  Eq.  (29),  twelve  terms  ( N  =  12)  gave  satisfactory 
results.  In  fact,  the  maximum  error  was  less  than  0.1%  between  the  applied  surface  load 
and  the  calculated  ones,  in  all  cases  computed.  The  maximum  error  occurred  very  near 
the  crack  tips  (x  -  0.99c).  Elsewhere,  along  the  crack  surface,  boundary  conditions  were 
satisfied  with  much  greater  accuracy.  Table  1  display  the  degree  of  accuracy  by  compar¬ 
ing  the  calculated  surface  load  with  the  applied  one  in  several  points  along  the  crack  sur¬ 
face. 
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The  error  in  boundary  tractions  near  the  crack  tips  is  the  result  of  the  approximation 
(32)  made  in  solving  the  dual  integral  equations  (25)  and  (26).  In  previous  work,  Eringen 
and  his  co-workers  have  shown  that  such  an  approximation  does  not  cause  appreciable 
errors  in  the  stress  field.  The  stress  field  calculation  uses  the  full  value  of  the  function 
K  (pq)  in  Eq.  (25). 

In  all  Figures  we  give  contributions  of  dislocations  and  dislocation-crack  interaction 
to  the  stress  field,  separately  from  the  total  stress  field.  The  total  stress  fields  acquire 
maxima  near  crack  tips;  just  outside  the  crack  tips.  Contrary  to  the  classical  elasticity 
solution  there  is  no  singularity  either  at  crack  tips  or  at  the  points  of  application  of  dislo¬ 
cations.  Consequently,  a  yield  criterion  may  be  set  up  by  equating  the  maximum  stress 
to  the  cohesive  stress  that  hold  bonds  together.  Such  a  yield  criterion  would  be  more 
meaningful  for  a  distribution  of  dislocations  along  the  slip  lines.  Such  an  investigation 
requires  extensive  analytical  and  computational  efforts,  extending  our  previous  study 
(13J  on  Mode  HI  crack  problems. 
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Table  1. 


Comparison  of  Applied  and  Calculated  Surface  Tractions. 


e  =  1 

a  = 

0-5  cm.,  c  = 

5  x  10-5  cm., 

(Fig-  3) 

10_  3  cm. 

a  =  67° 

e  =  10-5  cm.,  c  =  10~3  cm. 

a  =  10" 4  cm.,  a  =  75° 

(Fig.  6) 

x/c 

Applied 

Calculated 

Applied 

Calculated 

0.010 

-1.35830 

-1.35717 

-1.35238 

-1.35244 

0.248 

-1.36520 

-1.36486 

-1.35895 

-1.35896 

0.500 

-1.39340 

-1.39502 

-1.38573 

-1.38594 

0.752 

-1.49593 

-1.49412 

-1.47954 

-1.47953 

0.990 

-3.75938 

-3.76026 

-1.80860 

-1.80861 
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Fig.  1  Crack  and  four  dislocations 
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